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Section  1 
INTRODUCTION 

The  primary  purpose  oi  this  document  is  to  ))resent  a  generalized  method 
loi'  pi'odicting  the  blast  and  thermal  cnviionment  resulting  1  x-om  explosions  of 
tlie  three  most  common  liquid  propellant  combinations: 

•  Liquid  oxygen/RP-l  (LO^/RP-l) 

•  Liquid  oxygen/liquid  hydrogen  (LO^/LH^) 

•  Nitrogen  tetroxide/50%  liydrazine  -  50%  ue.symmetrical  dimethylhydrazine  - 
(N20^/50%  N2H^-50%  UDMH) 

The  blast  prediction  method  takes  into  account  the  fact  tliat  the  explosive 
potential  of  a  given  liquid  propellant  combination  in  accidental  failures  is 
not  a  unique  value,  but  depends  on  the  manner  in  which  the  propellants  are 
brought  together  during  the  failure  process  and  on  the  time  of  Ignition.  Spe¬ 
cifically,  the  method  provides  a  means  for  predicting  the  blast  environment  ^as 
a  function  of  the  following  system  characteristics: 

•  Tank  configuration 

•  Propellant  type 

•  Propellant  weight 

•  Nature  of  failure  mode 

•  Ignition  time 

The  thermal  environment  prediction  is  limited  to  the  characteristics  for 
each  of  the  propellant  combinations,  partly  because  of  the  more  limited  data 
available  and  partly  because  the  thermal  characteristics  are  less  dependent  on 
the  details  of  the  failure  conditions  than  are  the  blast  characteristics. 

It  should  be  emphasized  that  the  puipose  of  the  prediction  method  presen¬ 
ted  in  this  section  is  to  provide  a  means  for  predicting  the  explosion  environ- 
ment  as  a  function  of  the  various  controlling  parameters.  As  noted,  the 
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cryogGiiic  propellant  combinations  can  have  a  wiiJo  rniiKc  of  explosive  yield  val¬ 
ues,  depending  on  the  manner  in  which  the  propellants  come  into  contact  with 
each  other  and  mix  and  on  mixing  time  available  before  ignition.  Thus,  in  order 
to  use  the  prediction  method  for  any  given  vehicle,  it  is  necessary  to  eoiuluet  a 
detailed  failure  mode  analysis  to  establish  the  credible  failure  modes  and  the 
credible  ranges  in  values  lor  the  controli  ing  paiaiiietcrs . 

Included  in  this  volume  ai'e: 

Section  2.  A  brief  summary  of  the  general  scojie  and  the  types  of  results 
obtained  from  the  Project  PYRO  program,  (For  moi'e  details  the  reader  is 
referred  to  Volume  1,  the  Tlschnical  Documentary  Report  on  Pioject  PYRO.) 

Section  3.  The  blast  prediction  method  for  the  cryogenic  propellant  com¬ 
binations  , 

Section  4.  The  blast  piedictlon  method  for  the  hypergolic  propellant  com¬ 
bination  . 

Section  5.  Examples  demonstrating  the  use  of  the  cryogenic  blast  predic¬ 
tion  mjthod.  Whei'e  possible  the  cases  have  been  selected  to  correspond  to 
actual  full-scale  or  large-scale  failure  incidents  or  tests,  In  all  cases 
where  comparisons  are  made,  the  experimental  results  were  not  used  in  the 
basic  derivation  of.  the  prediction  method. 

Section  6  ,  Themal  p.redlction  method. 

Appendix  A.  A  discussion  of  liquid  propellant  explosion  pl'-’nomona  in  I’o- 
l&tlon  to  TNT. 

Appendix  B.  TNT  reference  overpre.ssure  and  impulse  versus  distance  curve.s. 
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Sod  ion  2 

IIRVILW  OF  PROJECT  I’VRO 

The  objective  of  Pi-ejccl  FYRO  wiis  to  develop  n  relinble  philosopliy  lor  pre- 
uidirij,  the  eieuiblc  damage  poLeiitinl  which  may  be  exj’er lenecd  from  the  neeiden- 
tnl  explosion  of  liquid  propellants  during  launch  or  lest  operations  of  military 
missiles  or  space  veliicles. 

Such  infojination  is  required  for  the  siting  of  static  test  stands  ami 
launch  facilities,  for  'iiifining  hazard  envelopes,  for  launch  operations,  etc. 

Tlie  PYRO  program  included  exi)erimental  determination  of  the  blast  and  ther¬ 
mal  environments  resulting  from  various  types  of  pi-opellant  mixtures  for  the 
throe  liquid  propellant  combinations;  N.^0^/5y%  -  50%  UDMH,  LO^/RP-1  ,  and 

LO2/LH2.  Propellant  weights  up  to  about  100,000  lb  were  used  for  the  cryogenic 
combinations  and  up  to  1,000  lb  for  the  hypergolic  combination. 

The  genoialized  test  conditions  used  wore  selected  to  simulate  the  impor¬ 
tant  classes  of  propellant  Interaction,  l.e.,  the  manner  in  which  the  two  pro¬ 
pellants  come  into  contact  wiiii  cacli  other  and  mix  during  an  accidental  failure. 
The  ways  in  which  the  proijollanls  can  mechanically  interact  witli  each  oilier  arc 
dependent  on  the  initial  conditions  of  the  propellants  at  the  start  of  the  in- 
ter.action  and  on  the  nature  of  the  boundary  conditions  wl,ich  eontiol  or  confine 
the  flow  of  propellants  during  the  spillage  and  mixing  process.  Tlie  two  major 
boundary  conditions  selected  for  testing  were  confinement  by  the  missile  and 
confinement  by  the  ground  surface. 

Tlie  conf Inemont-by- thc-missile  (CBM)  condition  is  intended  to  simulate  the 
general  case  where  failuie  occurs  in  the  intertank  bulkhead  and  all  jiropellant 
mixing  is  confined  within  the  tankage.  The  Initial  conditions  of  primary  con¬ 
cern  for  this  case  are  the  size  of  the  opening  in  the  intertank  bulkhead,  the 
length- to-diameter  ratio  of  the  tankage,  the  ullage  volume,  and  the  jiressure 
rise  to  cause  tank  rupture. 
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Tlig  conllnemciit-liy- tlic-trouncl-surfatc  (C1K5S)  Loiulltloii  .simulates  tin;  Lu.sg 
where  the  propellniUs  spill  out  of  the  tonknite  onti  mix  i)ii  the  r.rovnni  suifiice. 
Major  emphasis  in  the  program  has  been  planed  on  a  flat  pioiind  sarlnee,  al¬ 
though  a  limited  amount  of  data  was  obtained  for  other  conditions.  The  initial 
conditions  of  primary  concern  for  this  bouadaiy  condition  are  the  velocities, 
shapes,  and  relative  orientation  of  the  propellant  masses  at  the  start  of  their 
Interaction. 

Conceptually,  propellant  mixing  con  also  occui  without  confinement,  l.c., 
after  the  propellants  spill  out  of  the  tankage  but  befoi'e  they  rcacli  the  ground. 
Such  free-fall  mixing  was  not  included  in  the  program,  however,  becau-sc  of  llie 
small  amount  of  mixing  anticipated.  Unless  there  is  a  large  velocity  difference 
between  the  two  propellants  (which  is  unlikely  for  massive  fallui-es  neaj-  the 
gix>und),  there  are  no  significant  forces  holding  the  two  masses  together,  and 
even  a  small  pressure  generated  by  vaporization  or  reaction  nt  tlic  interface 
between  the  two  maoses  will  be  sufficient  to  separate  thrru  and  minimize  mixing. 

Although  the  generalized  test  conditions  used  i-esemble  some  actual  fi.ijlure 
modes,  the  intent  in  the  program  was  not  to  investigate  all  credible  eumblnaiions 
of  tankage  configuration,  failure  mode,  and  site  geometry;  there  is  an  almost 
infinite  number  of  such  combinations.  Rather,  it  was  reasoned  that  many  of 
these  combinations  would  lead  to  similar  propellant  interactions  and  that  study 
of  several  basic  propellant  Interaction  modes  would  be  sufficient  to  provide  a 
basis  for  evaluations  or  predictions  under  a  variety  of  failure  conditions. 

The  basic  blast  data  obtained  from  these  tests  were  peak  overpre.ssure  and 
positive-phase  impulse,  both  as  a  function  of  di.stance  from  the  propellant  ex¬ 
plosion.  Equivalent  explosive  weights  at  each  measurement  distance  were  deter¬ 
mined  separately  for  peak  overpressure  and  positive-phase  impulse,  using  stan¬ 
dard  TNT  surface  burst  reference  curves.  Characteristically  the  TNT  equivalent 
weights  computed  from  these  data  vary  both  as  a  tunction  of  the  shock  wave 
parameter  used  (peak  overpressure  or  positive-phase  impulse)  and  tlie  distance 
from  an  explosion.  At  long  distances,  however,  the  equivalent  weights  tend  to 
approach  an  equal  and  constant  value,  which  has  been  defined  as  the  terminal 
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cqui  vnlciil  wcluhl  (wlwti  t'xprcs.stHl  In  [louii'ts  of  ’I'NT)  oi  l(^'iiinin\l  yield  (when  ex- 
lucbsed  as  a  pereenl  ol  the  total  p lope  1  1  a ii t  welpht). 

The  yie’ld  values  based  on  peak  ovei'pressuie  generally  teiul  to  ineieiise  with 
Ineieaslni;  distance  I  roiii  the  explosion  until  the  terminal  yleltl  value  Is  renelied, 
Hie  yield  values  based  on  pus  i  1 1  ve-pliase  impulse  i;enerally  tend  to  decrease  with 
Inereaslnp  distance  until  the  terminal  yield  value  is  retielied. 

Thermal  data  obtained  if om  the  tests  included  total  heat  ilux,  nas  temi>era- 


lures  ,  aiKi  rudianl  he.,  t  ;  '  j 
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Section  3 

PHYO  HI, AST  I’liKDlCiTON  NttlTIIOI)  KOK  1,0., /Ill'- 1  AND  I,U./LII,, 

In  sctlll\^;  up  the  p  icii  1  c  I  ion  ihciIkhI,  it  nils  rocoi'.u  I  ■'cul  that  there  would 
likely  be  a  wide  eariety  oi  useis  (or  tiu’  method  and  that  tlie»’e  would  be  \ai  i- 
ou.s  level.s  oi  detail  deslvei,!  in  t  lie  output.  Aeeortli  nplj  ,  the  ha.sje  procedure 
will  deyiitned  to  piovide  tiie  mo.st  <letailed  intoniiation  possible  about  Ihi'  blast 
environment  within  the  current  state  ol  knowleilpe  repardinK  liquid  propellant 
e.\plosiVes,  but  in  a  lonnat  that  would  also  lieimll  le.ss  detailed  iniormation  to 
be  obtained  in  n  convenient  fashion. 

For  the  puiposes  of  predictint;  the  loaitinj;  and  response  of  objects  to 
bla.st  waves,  it  is  normally  .sufficient  to  specify  the  peak  overpressure  (!’)  and 
positive-phase  iminilsc  (I)  in  the  blast  wave. 

Tlius  the  basic  objective  ol  the  pretliction  method  was  to  provide  a  niea  is 
for  de  termini. ig  P  and  I  ns  a  lunctittn  of 

distance  -  U 

propellant  weipht  -  W 

characteristics  ol  mis.sile  or  space  vcdiiclc  .system  such  its  propellant  t.\  pe , 

tankage  (and  test  .site)  coni  1  gu  tit  t  ion ,  and  lailure  mode. 

The  basic  stejis  in  ttsing  the  prediction  method  are; 

1.  Conduct  engineering  analysis  of  specilic  sjstem  to  determine  ereilible 
failure  modes. 

2.  Dclennine  for  each  ereilible  failure  mode  the  general  seipieiice  of  events 
up  to  time  of  Ignition,  including  gross  space-time  history  of  propel¬ 
lants  . 

3.  Biiscd  on  the  itostulated  (tr.ipellant  space-time  history,  select  the  ;tp- 
[troitriate  generalized  PYKO  mixing  mode. 


In  the  remainder  oi  this  volume,  "peak-ovet  pressure"  and  ''iH>.si  t  i  ve-phase  im¬ 
pulse"  are  referred  to  ns  "overpressure"  and  "impulse." 
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tic, 

■J  ,  I,  ^  t  1111. 1  li_  c  .\|M  t  1 1  (I  \  11  I  lu  -  s  t  ‘  1  1 .1111, 1  ■■  "1  V  ;i  1  111.  )  'I  L'  111  1  I  "  1  1  1  lit;  I  11 1  .1  m- 

clLT'.  Im  ciiili  Mii.\liii;  tiioilf 

.  ttil<.r  id'I'liCs)  I'll  .Tiipi '  ii'i  1  :i  li  I'l  ■  ipt  1  1  .111  I  l>|’<  iiiid  inixiii).  iin)d'.-  n  1  t  li 
valucM  <il  olliL'i  I'll  i  .11111  1 1- 1  s  nipl  iil.ilain  iLimlniil  I'lld  iiiIul-.'^, 

0,  tnli!  I  (■  1 1.  IL  111  1.  TNI  p  I  hii  I  «-■— II  lid  1  nipvi  1  SI, -^1  1  s  I  ii  111  i.'  lUiMs  ()  lii.'ii  1  ii 

AppiliillX  b),  iislni;  t  q'l  1  vii  Kill  TN'I  '‘clslit  df  It  i  iii  I  lied  linii  Kinililili  >li  Id 
111  iiliiiiin  pitssiiit  (iiid  liiipiilst  lal'its  at  a  p.iitii  il  i  s ‘.iiiu  t  . 

7.  tiilti  .oliK's  lui  iippi  opi  i.'i  It  pi  iipt  1  i -'I  II 1  I  J  pt  wllli  (.tnltd  disUiiiLt  x.ii'it 
and  obtain  nppi  up  i  1  ii  It  to  i  itt  I  loii  l.actoi  loi  imiuilst.  (Kof  wtnltd  uls- 
taiitts  less  than  2  It'll)  stt  >|iia  1  i  I  >  i  I'p,  tuniiiit  ii  I  s  a  i  tiul  ol  Appiiidix 
A.  ) 

Tilt  pi  td  1  f  1 1  on  iTittliol  Is  piosfiiltd  In  a  scilt.s  ol  tbaits  (Ki|;s,  3-1  tliitninb 
3-r.).  HKiiit  3-1  blioMs  tilt  rtlalluii  bt'twtfii  I  ht  i;tiit  i  a  1  i  .'til  I’V  KO  nilMiit;  modts 
and  lypl'iiil  lull-.scnlt  lolKiit  liiodfs  to  assist  In  tin.'  stlttllni;  ol  t  lit  nppio- 
prla I f  nil X  1  ni;  mode  . 

Tilt  tipi.!  I  ions  lolntiiii;  t'Xpl  o.s  i  vt  jltld  to  ignlllnn  1  Imv  in  the  prcditllon 

nittlioil  wfiv  t.N  I  nipolultd  to  lull  si'ali'  by  tisiiij;  the  lollosliiK  posliilaltU  stnl- 

iiii;  rcliit  lonsli  i  p .  Ex|iloslvt  yield  (In  ptrctiit  o|  the  total  propellant  si'lrhi) 

is  indopfiulciit  ol  tilt  propellant  Mciplii,  providiiiK  llie  time  ol  Ipnltlon  Is 

1  /  3 

•sealed  b,\  the  cube  root  ol  itie  propellant  weight,  l.e,  ,  y  =  1  (l/W  '  ).  Tills 

.•■ealiiip  re  la  I  ioiislilp  was  po.stulnted  on  t  !•')  basis  ol  the  best  available  Inloniia- 
lioii  coiici'riiiiig  the  physical  plieiiomena  involved  in  the  propellant  mixing  and  ex¬ 
plosion  pi'ocessos  .  Tlie  Initial  check  ol  scaling  was  made  belwoen  tost  results 

1  roni  the  2UU  and  l,0du  ib  .scales,  and  the  I'esiilt.s  lor  all  cases,  c.xcept  those 

* 

lor  the  LO.^/RP-l  condition,  were  eonslstent  with  the  inistulalcd  scaling. 

(Thi'se  conijia ri sons  involved  200-  and  l,t)0t)-lb  tests.)  Later  a  limited  number 
ol  tests  weie  eoiidueled  in  the  range  Ironi  25,t)t>0-  to  1 00  ,Utl(J-lb ,  which  also 
showed  consistency  with  the  postulated  scaling.  These  larger  scale  comparisons 
arc  given  in  Section  5, 


*  Foi-  the  L0.)''Kr-l  CBM  case,  the  results  indicated  a  imiderate  ilecrcnso 
with  increasing  weight  up  to  aUnit  10,000  lb,  Abi'Ve  ihl.s  weight  the 
were  consistent  with  the  postulated  scaling. 


i  n  y  i  e  1  d 
resit  Its 
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mm  Fuu-scuE 
FAILURE  MODES 


PYRO  MIXINO  MOOES 


PREDICIION 

PROCEOURE 


BULKHEAD  RUPTURE 


Tailuro  occurk  in  the  intertank  bulkhead 


► 


9-id  oM  propellant  mixing  is  confined 
wi:nin  the  tankage. 


cam  •  CONFIN{Mf:N1  BY  MISSILF 
Cont'olling  ^oiametert; 

T  PROPfMAIMT  MIXING  TIME  -  msec 

lime  between  stort  of  mixing  ond  ignition, 

W  TOTAL  PROPELLANT  WEIGHT  -  lb 
t,  D  LENGTH-to-DIAMETER  RATIO  Of  PROPELLANTS 

0  '0  INTERIANK  BULKHEAD  OPENING  RATIO 
°  kotio  of  diameter  of  opening  in  intertank  bulkhead 

tc  lonk  diometer . 

V  TANKAGE  ULLAGE  VOLUME 
^  Percent  oi  totol  tankoge  volume.*'^ 

AP  TANK  RUPTURE  PRESSURE  DIFFERENTIAL  -  pii 
Burst  pressure  minus  initial  pressixe.U! 

^’*Not  needed  for  LC^/LHj. 


► 


LOj  'RP-1 

LOj/'lHp 


■■  See  Fig,  3'2 
-  See  Fig  >3 
(top) 


OVERPRESSURUATION-SEAM  RIP 
AND  low-altitude  FALLBACK 


► 


caos  «  confinement  6T  ground  surface 


► 


LOj  ^RP'l 
LOp/LHj 


See  Fig,  3-4 
See  Fig.  3-5 


Propellunts  spill  out  of  the  vehicle  and 
mixing  occurs  on  (he  ground  surface  . 
Impact  velocities  <  140  ft/sec. 


Cont'olling  Porofsseters : 


I  PkOPElLANf  MIXING  TIME  -  sec 

Tirrse  between  propellant  contact  on  ground  surface 
orsd  ignition. 

Vv  TOTAL  PROPEl  I  ANT  WEIGHT  •  lb 
V  PROPELLANT  VELOCITY  -  ft/sec 
FT  FAILURE  SUBTYPE 


•  Moisive  tankoge  rupture; 

propeHonts  reltosed  through  cssentiolly 
Full  cross  lection  of  tankage,  both 
tanks  opened  at  about  the  some  time. 


•  Some  os  I,  but  releose  of  top  propellant 
fr‘  m  tonkoge  is  considerably  delayed 
OAer  bottom. 

e  Propcilorits  released  rhrough  openings 
significantly  smoMer  than  tank 
cross  section. 


POWERED  IMPACT 

AND  HIGH -ALTITUDE  FALLBACK 


► 


HVI  •  HIGH-VELOCITY  IMPACT 


Impact  velocities  from  140  to  600Ft/*.*c, 
Ignition  -Cisumed  to  occur  ot  impact. 


Controlling  Porometers; 


V  IMPACT  VELOCITY  -  ff/»ec 
OS  NATURE  OF  GROUND  SURFACE 


Hord:  Essentially  rso  penetrotion  of  surfoce 
by  impacting  tankoge. 

(Example:  concrete  ond  rock.) 


Soft-  Essentially  complete  penetration  of  ssxfoce 
by  impoettng  tonkoge. 

(Example:  water.) 


► 


All 

Propellants  -  See  Fig.  3-3 
(bottom) 


Fig.  3-1.  Relation  Between  Typical  Full-scale  Failure  Modes 
and  P\F.O  Mixing  Modes 
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10,/RP-1:  CBM 


slepi 


DLTfKMiNE  P  A  R  A  M  E  T  t  R  V  A  L  UJ- S 


T  PROPELLANT  MIXING  TIME  - 

Tinu'  Ijctwcen  stQit  of  mixing  uiitj  ignitioti 

I  If  ufiknow?i  conbicJei  T  ) 
moK 

V.  TOTAL  PROPELLANT  VvElGHT  -  It. 

L  i)  LENGTH  to- DIAMETER  RATIO  OF  PROPELLANTS 

0  [>,  INTEKTANk  BULKHEAD  OPENING  RATIO 

of,. 

Ratio  of  dionictfi  of  opcnnig  In  InfeiiooR 
I'ulNicciJ  to  loi>L  dir.TicIci  . 

'If  opeiiitig  •Jlomett'i  oot  ciicului  diomclei 
of  CMcIo  lioviiig  son^e  opct'ing  oiccjl. 

TANKAGE  ULLAGE  VOLUME  -  peiceiit 
of  total  toiikogc  volume 

AP,  TANK  RUPTURE  PRESSURE  DIFFERENTIAL  -  pii 

(burst  pressure 
minus 

initial  piessuie) 


step  7 


D  L  T  C  R  m  i  N  E  PROCEDURE 


PROCEDURE 

W  -  lb 

D, 

TIME 

A 

> 10,000 

•'  0.45 

Unknown 

B 

^10,000 

<  0.45 

Known 

C 

■  10,000 

-  0  45 

Either 

c 

any  volue 

^0.45 

Eilhei 

procedure 

A 


step  2 


OBTAIN  TERMINAL  YIELD  (Y) 


Us  mg  L  D  value, 

obtain  from  curve  A,  Figure  A-). 

-  terminal  yield  at  tank  rupture 
for  -  10%,  =  85  psi. 


% 


Upper  Bound  90% 
Q  B  /''gdictlon  Interval 


•  With  Vy  and  "^P^  values,  determine  k  value  fri  ,T 
Procedure  A. 

With  k  end  L/D  values,  enter  Figure  B-1  and 
ollowoble  T/W  ^  volue .  Compare  with  T/V/ 

1  3 

less  than  5  msec/'lb  ^  , 


Figure  B - ) 


f 


'  Enter  I 
and  wl 
deterT 

TERMI 


If  W  <10,000  lb.  and  D  /D*<0.45, 

O  ' 

obtoin  yield  value  Y  from  generol  predlctic 
Procedure  A  or  B  as  appropriote. 


21 


terminal  yield  is  given  by:  Y(1  ■* 
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UEE 

URS  SYSTCMS 

cor»x>RAiJor^ 


step  4 


D  F  T  t  R  M 1  N  E 

PRESSURE  AND  IMPULSE  VALUES 


Determine  equivalent  TNEWTEV^^j) 

Went  Yx  w  to.  y^qx  w) 

; 

^  3 

Compute  scoloi-l  i^isfonce  (D  W 

i 


PRESSURE 
Defermltte  Py 
suifocc  burst  cuive 
In  Appendix  B 


IMPULSE 

Deterniine  from  TMi 

suiface  buist  curve 
In  Appendix  B 


Compute  I  b/  multiplying  Ijisil" 
by  Impulse  Correction  Factor 
(See  Table  1 ) 


TABLE  1 

SCALED  DISTANCE 
RANGE 


CORRECTION 

FACTOR 

1,3 

2,0 


LOj/LH^:  CBM 


1  DETEkM  NE  PARAMETER  VA'.UES 


T  PROPELLANT  MIXING  TIME  -  msec 

Time  between  start  of  mixing  and  ignition 
(If  unknown,  consider  ). 

W  TOTAL  PROPELLANT  WEIGHT  -  lb 
l/D  LENGTH-to  DIAMETER  RATIO  OF  PROPELLANTS 

D  /D  INTERTANK  BULKHEAD  OPENING  RATIO 
®  Ratio  of  diameter  of  opening  In  Intertank 

bulkheod  to  tank  diameter. 

(If  opening  diameter  not  circular, 
use  diameter  of  circle  having  same 
opening  area ) , 


N£  PROCEDURE 

PROCEDURE 

L/D 

%  /D, 

A 

>5.0 

All  values 

A 

»  1.8 

S  0.45 

B 

s  5.0 

s-  0.45 

HIGH-VELOCITY  IMPACT:  all  propellants 


stepj 


DETERMINE  PARAMETER  VALUES 


V  IMPACT  VELOCITY  -  ft/sec 

GS  NATURE  OF  GROUND  SURFACE 

Hard;  Essentially  no  penetration 
of  surfoce  by  impacting 
tankoge.  (Example 
concrete  and  rock). 

Soft:  Essentlolly  complete 

penetration  of  surface  by 
impacting  tankoge. 
(Example  -  wntei ). 


''"2  O.T 


AIN  TERMINAL  YIELD  (Y)  AND  UPPER  90%  PRE 

- j 

From  the  graphs  of  Terminal  Yield  vi  Impact  Velocity  for  L02/RP-1  or 
determine  TERMINAL  YIELD  (curves  A  and  C)  and  UPPER  90%  PREDIC 
for  appropriate  propellant  type.  I 


lOa/LHj 

150  : 


w 

o 

Ui 


100  : 


“  0  ■-  . 

0 


IMPACT  > 


AFni’L-TK-t)H-92 


UMS  SYSTEMS 
CORF  JiTATXJN 


ELO  (Y)  AND  UPPER  90%  PREDICTION  VALUES  (Y^g) 


I  I  Upper  Bouf^  90% 

^  ]  /  Prediction  Inferval  i 

»i  /r^ - - - : - i 

J Prediction  Equotion 

*  0  40  80  120  160 

T/W'/^  (miec/lb'/^) 


Dci'^rmine 
UPPER  90% 

PREDICTIOh.  VALUE 
from  curve  B . 

(There  i*  a  90%  probability 
that  the  yield  will  be  below 
this  value .) 


Est,  Upper 
"Bwnd  90% 


I  Estimated 
/I  Maximum 
h  Time 


^40  1 

Z  20 

I  1  ■ 

t;;  0  -  ■  - :  1 

0  10  20  30  40 


•  Determine 
UPPER 

PREDICTION  VALUE 
from  curve  B , 


TAV'/^  (msec/lb 


lELD  (Y)  AND  UPPER  90%  PREDICTION  VALUE  (Y,g) 

! 

Irield  V!  Impact  Veiocity  for  L02.^RP-1  end  LO^/LH^  , 

(curves  A  and  C)  and  UPPER  90%  PREDICTION  VALUE  (curves  0  arsd  0) 

h- 


Lst.  Upper 
90%  Bound 

"n 

c  I 

Prediction  Curve 


£st.  Upper 


Esf.  Upper 

LO2/LH2  90^c  Bound  Prediction  Cui ve 

D  C 

ISO  ,  ,  .  ,  .  .  Eit.  Upper 

S.  vv  '■f  90%  Bound. 

9  100  ,  ,  Prediction  Curve 


»  1  90%  Bound  ^  jg 
Prediction  Curve  5  P 


step  4 


DETERMINE 

PRESSURE  AND  IMPULSE  VALUES 


Determine  equivalent  TNT  ^ 

Wint  -  Y,W  (orYygXV/) 


Compute  Scaled  distance  (D/^  TNT^ 

i 

PRESSURE 

Determine  Pffsjj  from  TNT 
surface  burst  curve 
in  Appendix  B 


Determine  Ijfyjy  from  TNT 
Surface  burst  curve 
in  Appendix  B 

X 


Compute 
by  Impul 

1  by  multiplying  l^|.^.|. 
se  Correction  Fortor 
(See  Table  I ) 

SCALED  DISTANCE 

CORRECTION 

RANGE 

FACTOR 

t  s  1/3 

(ft/lb 

1 .4 

2.0 

0  200  400 

IMPACT  VELOCITY  (ft/sec) 


Fig.  3-3. 
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LO2/RP-I:  CBGS 


stepj 


DETERMINE  PAkAMLiER  VALUES 


T  PROPELLANT  MIXING  TIME  -  vec 

(Tinir  between  proncNonf  contoct  on  ground 
Surface  and  Igtiition}. 

-  fop  propcUont 
Tg  -  bottom  propellant 

If  unknown,  consider  T 

'  max 

W  TOTAL  PROPELLANT  V^EIGHT  -  lb 

V  PROPELLANT  IMPACT  VELOCITY  -  ft. 

FT  FAILURE  SUBTYPE 
I  Fallback 

Massive  tankage  rupture, 
propellants  relcosed  through  essentially 
full  cross  section  of  tankoge,  both  toriks 
opened  at  about  the  same  time. 

T  To  ,  V  velocity  of  interfoce 
between  ptopellanfs 


'B  ' 


II  Fallback  • 

Samp  os  I,  but  release  of  top  propellar^t 
from  tankage  Is  considerably  delayed 
over  bottom . 

III  Non-Mossive  Tank  Rupture 

Propellants  released  through  openitigs 
slgriificontly  srnaller  than  tonk  cross  section. 


step  2 


determine  procedure 


procedure 


c 

D 


TIME 

Known 

Known 

Known 

Unknown 


failure 

SUBTYPE 

1 

II 

III 

I,  II,  III 


OBTAIN  TERMINAL  YIELD  (Y)  AND  UPPER  90%  PRED 


WlO-j  ^  w, 

Compufe  from  following  equation;  Y  -  Yc  - 

o  W| 

where  .  .  . 

*  total  LO2  weight 
^  RP-1  “  of  overlapped 


KP-1- 


Wj  '  total  w( 

Y5  -  specific 
(right  o 
ofWLc 

•  Compere  yield  with  val'  e  given  in  Procedure  A 

(tjse  smollest  yield  value,  but  in  no  case  use  less  than  Procedud 

•  ESTIMATE  UPPER  90®^  PREDICTION  VALUE  by  increasing  yield  estimate  by 

using  Procedure  A  . 


1/3 

♦  Use  same  procedure  os  in  A,  but  use  in  ploce  of  W  for  computing  T/W  ' 

weight  of  both  propellants  on  the  ground  ot  time  of  ignition.  i 

TERMINAL  YIELD  VALUE  is  multiplied  by  Wf  (rcther  than  W)  in  Step  4  to! 
In  no  cose  should  on  equivalent  weight  in  lb  be  used  that  is  less  than  give| 


DROP  HEIGHT  (ft) 


•  With  V  value, 
detennine 
TERMINAL  YIELD 

from  curve  A. 

This  is  Iho  maximum 
Terminal  Yield, 
•  e.  highest  value  for  any 
ignition  time . 


120  I 
100 
80 
60 
40 
20 
0  ■ 


25 


100 


// 

'M 

c 

— 

— 

i 

\w 

!  V 

— 

— 

1*  1  ■ 

0  20  40  60  80  100 

IMPACT  VELOCITY  (ft/sec) 


•  VI 
d 
U 
f| 


/  i 


<1 


1/3 

c  In  piece  of  W  <*ipwting  T/NV  ^  -^vhere  is  the  estimated 
I  grourxJ  at  time  of  *91  ). 

tiplled  by  Wf  (rather  ^  V)  in  Step  4  to  get  equivalent  weight  in  lb  . 

^eight  In  lb  be  used  tncT  1$  less  than  given  In  Procedure  A  with  a  LAA/  '' 


DROP  Hill,  ;T  (ft) 


?/  / 

100  i  ?  jV  / 

I  J~s/ 

T  jt/  /'^ 

i^/A 

i  ■-//// 


With  V  volue, 
determine 

UPPER  90%  PREDICTION  VALUE 
from  curve  C . 


0  20  40  60  80  100 

IMPACT  VELOCITY  (ft/sec) 
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UHS  SY5>TEMS 
CORt’URATlON 


StfiDil  i^etermine 

PRESSURE  AND  IMPULSE  VALUES 


Deternune  equivalent  TNT  ^ 


Compute  s  alcd  distance  (D  \V  '  ^TNT^ 

i 

PRESSURE 

Determine  P,,,-  from  TNT 
TNT 

suifoce  buisl  cuivc 

in  Appendix  B  _  » 


Determine  1-^,,  fnxn  TNT 
TNT 

surface  buist  curve 
In  Appendix  B 

4- 

Compute  I  by  multiplyifig 
by  Impulse  Coirectlon  Foctoi 
(.See  Table  1 ) 


TABLE  1 

SCALED  distance 
RANGE 

(fLTb’/ 

CORRECTION 

factor 

1.3 

2.0 

TABLE  2 

LO2  ^RP-I 
Rotio 

Ys  ("0) 

1 .5 

93 

2.0 

113 

2.5 

126 

3.0 

132 

3.5 

123 

4.0 

115 

4.5 

T05 

5.0 

96 

6.0 

83 

7.0 

70 

8.0 

59 

10.0 

43 

12,0 

31 

14.0 

22 

16.0 

16 

18.0 

12 

v 


ims  G52-35 


LO,/lK,:  CBGS 


stepi 


DETERMINE  PARAMETER  VALUES 


T  PROPELLANT  MIXING  TIME  -  »ec 

(Time  ber^een  propellanf  contaci  on  ground 
»urfac«  and  ignition). 


Tg  -  bottom  prc^ellant 

If  unknown,  consider  T 

mcK 

W  TOTAL  PROPELLANT  WEIGHT  -  lb 
V  PROPELLANT  VELOCITV  -  ft/sec 
FT  FAILURE  SUBTYPE 

I  Fallback  - 

Massive  tankage  rupture, 
propellonts  released  through  essentially 
full  cross  section  of  tankage,  both  tanks 
opened  ot  about  the  some  time. 

T  =  Tg  ;  V  -  velocity  of  interface 
between  propellonts 

II  Fallback 

Some  os  I,  but  release  of  top  propellant 
from  tankoge  is  considerably  deloyed 
over  bottom . 

Ill  Non-Massive  Tank  Rupture 

Propellants  released  through  openirsgs 
significantly  smaller  than  tank  cross  section. 


ERMINE  PROCEDURE 


PROCEDURE 

A 

B 

C 


TIME 

Known 

Known 

Unknown 


FAILURE 

SUBTYPE 

I,  It 
III 


procedure 

h 


steps 

OBTAIN  TERMINAL  YIELD  (Y)  AND  UPPER  90%  PREQ 


*  Co-npute  T/^  from  input  data 
and  dttefmine  TERMINAL  YIELD. 


•  With 

detftmine  UPPII 


^ - -7 

vtuH  r  1  EH.  t  -1 II,*'' 

»<  '  tl 

f/y 

y - 1 

T  tr  JF 

4  V  4.  , 

i  M 

_ i 

'  1/3  ’ 

*  At  T,^  •  0;  Y  -  18%. 

0.2 


0.3 


0.4 


T/W 


e  Use  same  procedure  as  in  A,  but  use  in  piace  or  vy  ror  comp' 
weight  of  both  propellants  on  the  ground  ot  time  of  ignition 


•  lace  of  W  for  computing 


TERMINAL  YIELD  VALUE  is  multiplied  by  (rather  than  W)  in 
In  no  cose  should  on  equivalent  weight  in  lb  be  used  that  is  less  ' 
with  a  T/W’'^^  =  0. 


•  With  V  voiue, 
determine 
TERMINAL  YIELD 
from  curve  A. 

This  is  the  moKimum  Terminal  Yield, 
i  .e.  highest  value 
for  ony  ignition  time. 


IMPACT  VELOCITY  (ft/sec) 


Fig.  3-5. 
3-11 
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Section  4 

PYRO  BLAST  PREDICTION  METHOD  FOR  N,.0  /50%  N.  -  50%  UDMH 

2  4  2  4 

The  general  approach  used  In  the  prediction  method  for  NjjO^/50%  -  50% 

UDMH  Is  similar  to  lliat  lor  the  cryogenic  propellant  combinations.  The  pixjce- 
dures ,  however,  are  considerably  simpler  because  time  of  ignition  is  not  a  con¬ 
trolling  parameter  and  because  the  generally  low  magnitudes  of  the  yields  made 
it  possible  to  eliminate  much  of  the  parameter  variations  needed  for  the  cryo¬ 
genic  propellent  combinations. 

The  basic  steps  in  the  prediction  procedure  are: 

1.  Identify  the  PYRO  generalized  configuration  from  the  list  given  below. 

2.  Follow  specified  procedures  to  determine  terminal  yield  and  pressure 
and  Impulse  correction  factors. 


GENERALIZED  CONFIGURATION 

1.  Static  Test  Stand 

a.  Assumes  tank  and  tank  support  structure  are  strong  enough  to  elim¬ 
inate  failover  case  and  no  large  explosive  donor  (i.e.,  another 
stage)  Is  present. 

b.  Applicable  PYRO  mixing  modes  are:  diaphragm  rupture  (confinement 
by  the  missile);  spill;  small  explosive  donor  (i.e.,  sO.5%  of  the 
propellant  weight),  and  tower  drop. 

2.  Launch  Pad 

a.  Pre-launcl» 

Applicable  PYRO  mixing  modes  are:  diaphragm  rupture,  spill,  small 
explosive  donor,  and  lower  drop. 

b .  Launch 

Applicable  PYRO  mixing  modes  include  those  indicated  for  pre-launch 
as  well  as  large  explosive  donor  (I.e,,  another  stage)  and  low- 
velocity  impact  (~140  ft/sec).  Docs  not  include  high-velocity  im¬ 
pact  (either  from  hlgh-altitude  fallback  or  powered  Impact). 
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3.  Pust-Launch 

a,  In-Fliglit 

No  data  wero  obtained  lor  this  case  duriiig  tlic  I’YRO  ijroKi'iiiri;  liow- 
ever,  In  the  absence  of  a  spill  surface,  there  seems  to  be  no  pos¬ 
sible  way  for  the  failure  mode  to  be  more  severe  than  the  launch 
case.  It  is  recommended  that  this  value  be  used. 

b.  Ground  Impact 

Appllcoble  PYRO  mixing  mode  is  hlgh-volocl ty  impact. 

■nible  4-1  gives  the  estimated  upper  limit  for  the  terminal  yield  for  eacli 
of  the  PYRO  generalized  configurations.  For  all  cases  but  the  hlgli-velocl  ty  im¬ 
pact  condition,  a  single  ternilnal  yield  value  could  be  u.sod.  For  this  case  the 
voluc  shown  In  Table  4-1  is  for  the  upper  limit  of  velocity  (lioo  fl/soc)  and  a 
curve  of  terminal  yield  vs  velocity  Is  given  In  Fig.  4-1  for  use  with  lower  im¬ 
pact  velocities. 

With  the  terminal  yield  for  n  specified  failure  mode  establisliod  by  use  ol 
Table  4-1  or  Fig  4-1,  the  next  step  In  the  prediction  procedure  Is  to  detcnnine 
the  pressure  and  Impulse  corjxjctlon  factors.  These  arc  obtained  from  Fig.  4-2 
by  entering  in  the  pixiper  scaled  distance.  Tlio  propellant  overpressure  and  im¬ 
pulse  values  are  obtained  by  multiplying  the  overpressure  and  impulse  values 
obtained  from  the  surface  burst  TNT  curves  shown  in  Appendix  B  by  these  factoi'S, 


4-2 


URS  002-35 


AFRPL-TK-68-93 


Table  4-1 

TERMINAL  VIELD  ESTIMATES  FOR 
SELECTED  FACILITY  AND  APPLICATION  MODES 


FACILITY  AND 

ESTIMATED  UPPER  LIMIT 

APPLICATION  MODE 

TO  TERMINAL  YIELD  (%) 

STATIC  TEST  STAND* 

2 

LAUNCH  f>AD 

* 

#  Pre-Launcti 

2 

•  Launch 

5 

POS'P- LAUNCH 

•  In-Flight* 

5 

•  High  Velocity  Impact 

Hard  Surface 

25 

Soft  Surface 

60 

If  no  small  explosive  donor  (i.e., 
present  a  yield  of  1.5  percent  can 

If  no  largo  explosive  donor  (i.e., 
3  percent  can  be  used. 


<1%  of  total  propellant  weight)  is 
be  used  . 

another  stage)  is  present  a  yield  of 


4-3 


TERMINAL  YIELD  (%) 


URS  652-35 


AFRPL,-TR-68-92 


IMPACT  VELOCITY  (ft/s«c) 


Fig.  4-1.  Terminal  Yield  vs  Impact  Velocity  for 
Kypergollc  High-Velocity  Impact 
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PEAK  OVERPRESSURE 


4  5  6  8  10  15  20  30  40  60  80 


SCALED  DISTANCE  (ft/lb  V3) 


POSITIVE-PHASE  IMPULSE 


4  5  6  8  10  15  20  30  40  60  80 


SCALED  DISTANCE  (ft/lb’/^) 


Fig.  4-2.  Ratio  of  Upper  Bound  of  Propellant  Peak  Overpiessui’e  and 
.Positive  Pha.se  Impulse  Data  to  Standard  TOT  Curve 
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Section  5 

examples  of  use  and  experimental  verification 

OF  PREDICTION  METHOD 


Examples  illustrating  the  use  of  the  prediction  method  are  given  in  this 
section.  The  cases  covered  are  listed  below.  Where  possible  the  cases  have 
been  selected  to  corrcsjjond  to  actual  full-scale  or  large-scale  failure  inci¬ 
dents  or  tests.  In  all  cases  where  comparisons  are  made,  the  experimental  re¬ 
sults  were  not  used  in  the  basic  derivation  of  the  prediction  method. 


EXAMPLE 

NO. 

PROPELLANT 

COMBINATION 

FAILURE  MODE 

COMPARISON 

CASE 

TYPE  OF 
COMPARISON 

1 

LO^/RP-l 

Bulkhead  Rupture  -  Ignition 
Time  Unknown 

ATLAS  9C 

Full-Scale 

Accident 

2 

LO^/RP-1 

Bulkhead  Rupture  -  Ignition 
Time  Known 

AFRPL  TITAN  I 
Test  (No. 

301) 

Full-Scale 

Test 

3 

LO/LH^ 

Bulkhead  Rupture  -  Ignition 
Time  Known 

AFRPL  SATURN 
S-IV  Test 
(No.  62) 

Full-Scale 

Test 

A 

L0.,/RP-1 

Fallback  -  Ignition  Time 

Known 

AFRPL  25,00C- 
Ib  Ibst 
(No.  285) 

Large-Scale 

Test 

5 

LO^/RP-l 

Fallback  -  Ignition  Time 
Unknown 

None 

None 

6 

LO/LH^ 

Fallback  -  Ignition  Known 

AFRPL  25,000- 
Tfest  (No. 

228C) 

Large-Scale 
Tes  t 

7 

LO2/RP-I 

LO^/LHa 

Fallback  -  Ignition  Time 

Known 

AFRPL  1  ,200- 
lb  Test 
(No.  295) 

Small-Scale 

Test 
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1.  BULKHEAD  RUPTltRE  OF  LO2/RP-I  VEHICLE  ON  LAUNCH  PAD  OR  STATIC  TEST  STAND  - 
IGNITION  TIME  UNKNOWN  -  CORRESPONDS  TO  PYRO  CBM  CASE 

♦ 

A:isumed  Conditions 


W  =  240,000  Ib  (weight  of  propellants  when  vehicle  has  normal  full 
load  of  propellants) 

L/D  =  5 

D  /D^  £  0.45 
o  t 

V  =  40%  at  time  of  failure 
u 

AP  =  30  psi 
r 

Prediction  Procedure  (Refer  to  Fig.  3-2) 

1.  In  accordance  with  step  2,  Procedure  A  is  selected. 

2.  Following  step  3,  a  Yg  of  10%  Is  selected  from  curve  A  of  Fig,  A-1 
with  tne  given  L/D  value  of  5.  Also  a  k  of  0.9  is  obtained  from 
Fig.  A-2  using  APj.  =  30  psi  and  a  Vy  =  40%.  The  predicted  termin¬ 
al  yield  value  of  9%  is  then  obtained  by  multiplying  0.9  times  10% 

3.  Ck>ntlnuing  with  step  3,  the  upper  90%  prediiction  bound  is  14% 

(0.9  X  15%). 

4.  Following  step  4,  an  impulse  correction  factor  of  1.3  is  obtained 

for  scaled  distances  of  greater  than  3  ft/lb^*^^  and  a  value  of  2 
for  scaled  distances  less  than  3  ft/lb^^^.  With  the  derived  ex¬ 
plosive  yield  of  9%  or  22,000  lb  of  TNT  (for  the  240,000  lb  of 
propellant)  this  means  that  for  actual  distances  greater  than  3  x 
22,000f/3  :::  84  the  impulse  values  obtained  from  standard  INT 

impulse— distance  curves  should  be  multiplied  by  1.3.  Similarly 
for  distances  less  than  84  ft,  the  impulse  value  should  be  multi¬ 
plied  by  2. 


W  is  total  propellant  weight;  L/D  is  length- to-diame ter  ratio  of  propellants; 
Dy/Dt  is  ratio  of  the  diameter  of  the  opening  in  the  bulkhead  to  the  tank 
diameter;  V  is  tankage  ullage  volume;  AP  is  tank  rupture  pressure  differ¬ 
ential  . 
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Full-Scale  Comparison 

The  case  selected  is  believed  to  correspond  to  the  Atlas  9C  failure. 
The  actual  ullage  volume  at  ignition  in  this  case  was  unknown,  but  propel¬ 
lants  were  being  pumped  out  of  the  vehicle  for  30  to  60  sec  prior  to  igni¬ 
tion,  so  that  the  ullage  volume  was  clearly  very  much  larger  than  the 
normal  5%,  The  measured  yield  from  the  failure  was  10-12%,  which  compares 
quite  favorably  with  the  predicted  value  of  9%  and  the  upper  bound  of  14%. 


2.  BULKHEAD  RUPTURE  OF  LO2/RP-I  VEHICLE  ON  LAUNCH  PAD  OR  STATIC  TEST  STAND  - 
IGNITION  TIME  KNOWN  -  CORRESPONDS  TO  PYRO  CBM  CASE 

Assumed  Conditions 

T  =  840  msec 

W  =  170,000  lb  (weight  of  propellants  when  vehicle  has  normal  full 
load  of  propellants) 

L/D  =  4 

D  /D,  S  0.45 
o  t 

V  =  50%  (corresponds  to  a  propellant  loading  of  about  55%  or  94,000 
lb  of  propellant 

AP  =  35  psi 

r 

Prediction  Procedure  (Refer  to  Fig.  3-2) 

1.  In  accordance  with  step  2,  procedure  B  is  selected. 

2.  Following  step  3,  a  k  value  of  0.9  is  obtained  from  Fig.  A-2  by 
entering  with  a  Vy  =  50%  and  a  APj.  =  35  psi.  Next  a  maximum  al¬ 
lowable  value  of  T/W^'^2  =  27  msec/lb^^^  is  obtained  from  Fig.  B-1 
by  using  the  k  value  of  0.9  and  an  L/D  =  4.  Since  this  is  ^arger 
than  the  value  computed  from  the  input  data,  840/(170,000)^'^  = 

15,  the  input  value  (being  the  smaller  of  the  two)  is  selected  and 
entered  into  Fig.  B-1,  with  an  L/D  -  4  to  obtain  a  terminal  yield 
=  6%. 

3.  Again  using  a  T/W^'^''  =  15  and  an  L/D  =  4,  the  upper  90%  predic¬ 
tion  value  =  8%  is  obtained  from  the  lower  part  of  Fig.  B-2 . 
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Full-Scale  Comparison 


The  case  selected  corresponds  to  the  Titan  I  test  conducted  under  the 
FYRO  program  (Thst  Number  301).  The  measured  yield  was  4%  based  on  the 
94,000  lb  of  propellant  used  or  2%  on  a  fully  loaded  vehicle.  Tlie  predic¬ 
ted  value  of  6%  would  .le  expected  to  be  conservative  (‘.e.,  somewhat  high) 
for  this  case  because  of  the  small  s  ze  of  the  opening  area  in  the  inter¬ 
tank  bulkhead,  which  amounted  to  only  about  2%  of  the  tank  cross-sectional 
area.  The  prediction  method  was  based  on  results  with  diaphragm  openings 
equal  to  or  greater  than  20%  of  the  tank  cross-sectional  area  end  was  ex¬ 
pected  CO  be  conservative  for  cases  having  much  smaller  openings. 

3.  BULKHEAD  RUPTORE  OF  LO2/LH2  VEHICLE  ON  LAUNCH  PAD  OR  STATIC  TEST  STAND  - 
IGNITION  TIME  KNOWN  -  CORRESPONDS  TO  WRO  CBM  CASE 

Assumed  Condition 

.  183  j 

W  =  91,000  lb 

L/D  =1.8 

D  /D^  ^  0.45 
o  t 

Prediction  Procedure  (Refer  to  Fig.  3-3) 

1.  In  accordance  with  step  2,  procedure  A  is  selected. 

2.  Following  step  3,  a  T/W^'^^  =  4  msec/lb^'^  is  computed  (183/91,000^^^) 
and  a  terminal  yield  of  6%  is  obtained  by  using  this  value  with 
curve  A. 

3.  An  upper  90%  prediction  bound  of  12%  is  obtained  from  curve  B  in 
the  same  manner. 

Full-Scale  Comparison 


This  case  corresponds  to  the  S-IV  test  conducted  as  part  of  Project 
PYRO  (Test  No.  62).  The  measured  yield  was  5%,  which  compares  favorably 
with  the  predicted  value  of  6%. 
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4.  FALLBACK  OF  LO2/HP-I  V1';HICLE  ON  LAUNCH  PAD  -  IGNITION  TIME  KNOWN  -  CORRE- 

sponDl-  to  pyro  cbgs  case 

4( 

Assumed  Condi lion 
T  =  0.465  .sec 
W  =  25,000  Ib 
V  =  44  It/sec 
FT  =  I 


Prediction  Procedure  (Refer  to  Fig.  3-4) 


1.  In  nccordance  witli  .stej)  2,  procedure  A  is  selected. 


2. 


Following 

25,000^/^ 


sloj)  3 
)  and  a 


ue  and  a  V  =  44. 


a  T/W^'^^  =  0.016  sec/lb^'^^  is  computed  (0.465/ 
terminal  yield  of  36%  is  obtained  from  this  val- 


3.  Continuing  willi  stej)  3,  an  upper  90%  prediction  bound  of  46%  Is 
obtained  in  a  similar  fashion. 


Lni'ge-Sca  le  Compni'ison 


Tills  I'.'ise  cori-espotu.s  to  a  25,000-lb  lest  (No.  285)  conducted  under 
the  PYRt)  pro  Irani.  The  measured  yieid  was  37%,  which  compares  very  well 
with  the  |iri'ili  led  value  ot  36%  and  upper  bound  of  46%. 


5.  FALLBACK  OF  L0,,/llP-l  VEHICLE  ON  LAUNCH  PAD  -  IGNITION  TIME  UNKNOWN  -  CORRE¬ 
SPONDS  TO  pyho“cik;.s  case 

Assumed  Condition 
W  =  25,000  lb 
V  =  44  ft/sec 
FT  =  I 


V  is  impact  velocity;  FT  is  failure  .subtype. 
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PrecHctl('n  Pi'occduro  (Refor  to  Fl^;.  3-4) 

1.  In  accoi'Cinnce  wl\.Ii  step  2|  procedui’e  D  is  stfloctt^iJ . 

2.  Followlntr  step  3,  a  terminal  yield  value  of  97%  is  obtained  by 
Uf  trij,.  a  V  --  4*1  ft/sec. 

3.  aint'nulng  .  tth  step  3,  an  upper  90%  prediction  bound  of  120%  is 
obtained  in  a  sliniiar  fashion. 


6,  FALLBACK  OF  LO2/LH2  VEHICLE  ON  LAUNCH  PAD  -  IGNITION  TIME  KNOWN  -  CORRE¬ 
SPONDS  TO  FYRO  CBGS  CASE 

Assumed  Conditions 
T  =  0.365  sec 
W  =  25,000  lb 
V  =  44  ft/sec 
FT  =  I 

Prediction  Procedure  (Refer  to  Fig.  3-5) 

1.  In  accordance  with  step  2,  procedure  A  Is  selected. 

2.  Following  step  3,  a  =  0.012  sec/lb^^^  is  computed  (0.365/ 

25,OOo1/3)  and  a  terminal  yield  of  22%  is  obtained  from  this  val¬ 
ue  and  a  velocity  of  44. 

3.  Continuing  with  step  3,  an  upper  90%  prediction  bound  of  42%  is 
obtained  in  a  similar  fashion. 


Large-Scale  Comparison 


This  case  corresponds  to  a  25,000-lb  test  (No.  288C)  conducted  under 
the  PYRO  program.  The  measured  yield  was  13%,  somewhat  under  the  predic¬ 
ted  value  of  22%  but  well  within  the  uncertainty  limits  for  this  case,  as 
evidenced  by  the  difference  between  the  expected  value  of  22%  and  the  up¬ 
per  90%  bound  of  42%. 
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7.  fallback  of  two-stage  vehicle  (UPPER  STAGE  LO2/LH2  -  LOWER  STAGE  LO2/RP-I) 
-  IGNITION  TIME  KNOWN  -  CORRESPONDS  TO  PYRO  CBGS  CASE 

Assumed  Conditions 

W  -  LO  /RP-1  stage  =  1,000,000  lb 
W  -  LO^/LH^  stage  =  200,000  lb 

T  =  5.4  sec 
V  =  44  ft/sec 
FT  =  I 


Prediction  Procedure  -  LO2/RP-I  Stage  (Refer  to  Fig.  3-4) 


1.  In  accordance  with  step  2,  procedure  A  is  selected. 

1/3  1/3 

2.  Following  step  3,  a  T/W  =  0.054  is  computed  (5.4/1,000,000  ), 

and  a  terminal  yield  of  93%  is  obtained  from  this  value  and  a  V 
=  44.  This  gives  an  equivalent  explosive  weight  of  930,000  lb 
<93%  X  1,000,000  lb) . 

Prediction  Procedure  -  LO2/LH2  Stage  (Refer  to  Fig.  3-5) 


1.  In  accordance  w'lth  step  2,  procedure  A  is  selected. 

1/3  1  '3 

2.  Following  step  3,  a  T/W  =  0.092  is  computed  (5.4/200,000  ), 

and  a  terminal  yield  of  37%  is  obtained  using  this  value  and  a  V 
=  44.  This  gives  an  equivalent  explosive  weight  of  74,000  lb 
(37%  X  200,000  lb). 


Combined  Yield  Predictions 


Sum  of  equivalent  explosive  weights  is  1,004,000  lb  (930,000  +  74,000) 
giving  a  terminal  yield  for  the  combination  of  84%  (100  x  1,004,000/ 
1,200,000). 

Experimental  Comparison 

No  full-scale  or  large-scale  data  are  available  for  comparison  v.ith 
this  case;  however,  one  1,200-lb  test  (No.  295)  was  conducted  as  part  of 
the  PYRO  program  with  conditions  properly  scaled  to  match  these  conditions 
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(weights  scaled  by  factor  of  1,000,  times  scaled  by  factor  of  10). 
measured  combined  yield  for  this  test  was  70%,  which  compares  well 
the  predicted  yield  of  84%. 


The 

with 
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Section  6 

HEAT  TRANSFER  HAZARD 

Tlic  flrobnll  generated  by  the  explosion  of  propellant  mixtures  can  consti¬ 
tute  a  hazard  primarily  through  heat  transfer  to  an  object  or  structure  Immersed 
In  it.  This  section  contains  a  description  of  the  expected  dimensions  and 
duration  of  the  associated  fireball  and  of  the  heat  flux  density  with  time  with¬ 
in  it,  each  as  a  function  of  tlie  quantity  of  pi-opellants  involved.  Remarks  in- 
dlcatinK  basis  and  limitations  of  each  prediction  are  also  included. 

Tlie  dimensions  of  the  fireball  depend  generally  on  the  quantity  of  propel- 
lant.s.  An  empirically  derived  expression  relating  the  fireball  dimension  in 
terms  of  an  equivalent  diameter  D  in  feet  to  the  total  propellant  (fuel  and  oxi¬ 
dant)  weight  W  in  pounds  for  the  propellant  combinations  of  L0_/RP-1 ,  LO„/LH  , 
RP-1/LH2/L0^,  and  N^H^  -  50%  UI»IH  is  given  by 

D  9.56  (6.1) 

where  the  estimated  standard  error  in  the  diameter  is  30%,  Equation  (6.1)  does 
not  always  provide  an  accurate  indication  of  the  maximum  dimenslon(s)  of  the 
fireball,  since  "in  those  instances  where  the  fireballs  were  markedly  asymmet- 
rical ,  attempts  were  made  to  estimate  equivalent  spherical  diameters.  An 

.indication  of  the  departure  from  the  diameter  given  in  Eq.  (6.1)  of  the  maximum 
dimensions  that  can  occur  is  provided,  for  instance,  by  the  Titan  test,  which 
involved  approximately  100,000  lb  of  LO^/RP-l .  The  maximum  horizontal  fireball 


Infonnation  permitting  the  evaluation  of  thermal  hazards  external  to  the 
fireball  through  radiant  energy  transfer  are  given  in  Section  6,  Volume  1, 

Equation  (6.1),  along  with  Eq.  (6.2)  below,  have  been  extracted  from  J.  B. 
Gayle  and  J.  W  Bransford,  Size  and  Duration  of  Fireballs  from  Propellant 
Explosions ,  NASA  TM  X-5312,  August,  1965. 
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dimension  from  this  test  was  estimated  to  be  from  800  to  1000  ft,  while  Eq  (6.1) 
indicates  diameters  of  approximately  400  ft. 

The  fireball  duration  t  in  seconds,  that  is,  the  time  over  which  fireball 
temperatures  persist  at  hazardous  levels  (excluding  residual  fires  of  unburned 
propellants,  which  tend  to  collect  in  ground  surface  depressions  or  structural 
confinements)  is  given  by. 


T  =  0.196  (6.2) 

where  the  standard  error  in  the  duration  is  84%. 

Curves  from  wh.ich  the  heat  flux  density  with  time  within  the  fireball  can 
be  obtained  for  a  given  propellant  weight  are  given  on  Figs.  6-1  and  6-2  for 
the  lAi^/RP-1  and  LO^/hH^  propellant  combinations,  respectively.  The  time 
in  these  figures  is  given  in  seconds  by, 

1/3 

T  =  C  W  '  (6.3) 

o 

for  a  total  propellant  weight  W  in  pounds,  with  a  value  of  C  of  0,113  for  LO^/ 
RP-1  (Fig.  6-1)  and  of  0.077  for  LO  /LH  (Fig.  6-2).  TVvo  curves  are  presented 

u  A 

in  each  figure.  One  is  the  "bounding  curve,"  which  is  an  estimate  of  the  upper 
bound  of  the  heat  flux  density  and  is  primarily  based  on  the  analysis  of  heat 
flux  density  data  that  were  obtained  from  eleven  25,000-lb  propellant  tests, 
five  of  LO  /RP-1 ,  and  six  of  LO  /LH  .  The  remaining  curve,  designated  the 

A  A  A 

"recommended  curve,"  is  superimposed  on  the  bounding  curve  until  a  time  t  — 

o 

given  by  Eq.  (6.3)  —  where  it  abruptly  decreases  to  zero.  The  recommended  curves 
are  also  based  primarily  on  analysis  of  the  data  from  the  eleven  25,000-lb  tests 


Data  from  which  the  heat  flux  density  may  be  evaluated  for  the  N204/50% 

N2H4  -  50%  UDMH  propellant  combination  are  extremely  limited.  Examination 
of  these  data  suggests  that  the  heat  flux  density  is  somewhat  less  in  magni¬ 
tude  than  the  bounding  curves  given  for  LO2/RP-I  and  LO2/LH2  in  Figs.  6-1 
and  6-2,  but  that  the  heating  durations  are  perhaps  somewhat  larger. 
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mentioned  above,  and  implicitly  contain  the  constraint  that  the  probability  of 
exceeding  the  cumulative  heal  flux  density  associated  with  the  recommended 
curves  (the  time  integration  of  the  heat  flux  density  from  time  equal  zero  to 
T^)  is  The  variation  of  the  heating  pulse  with  propellant  weight,  that  is, 

the  scaling  implicitly  contained  in  Figs.  6-1  and  6-2  and  Eq .  (6.3),  assumes, 
first,  that  the  duration  of  the  iieating  pulse  will  Increase  with  the  cube  root 
of  propellant  weight,  as  implied  by  the  empirical  relation  Eq.  (6.2)  and,  sec¬ 
ond,  that  the  heat  flux  density  at  a  scaled  time,  using  this  cube  root  time 
scaling,  will  be  invariant  with  variation  in  propellant  weight.  The  second 
statement  is  based  on  the  invariance  of  fireball  temperatures  (measured)  from 
scale  to  scale. 

No  account  has  been  made  in  the  bounding  or  recommended  curves  for  the 
emission  of  radiant  energy  from  the  surface  of  an  immersed  object,  and  this 
emission  can  substantially  reduce  the  transfer  rates  from  those  given  in  the 
curves  as  the  surface  temperature  of  the  object  becomes  a  significant  fraction 
of  the  fireball  temperature,  the  latter  being  typically  of  the  order  of  2300°K, 

A  reduction  occurs  similarly  for  the  convective  component  of  transfer.  Any 
corresponding  modifications  of  heat  transfer  from  the  curves,  however,  depend 
on  the  details  of  the  application  and  are  not  considered  here. 

Several  other  qualifications  of  the  bounding  and  recommended  curves  should 
be  noted.  First,  the  heat  flux  density  measurements  upon  which  the  curves  are 
primarily  based  were  obtained  from  instruments  that  were  fixed  in  space;  thus, 
a  modified  heat  flux  density  may  be  appropriate  for  objects  which,  for  example, 
become  prematurely  ejected  from  the  fireball  (due,  for  instance,  to  blast  wave 
forces) .  For  many  circumstances,  the  modification  would  be  a  reduction  of  the 
total  heat  transfer,  first,  due  to  the  tendency  to  reduce  the  time  that  an 
object  is  immersed,  and  second,  due  to  a  reduction  in  the  convective  heat 
transfer  component,  since  the  motion  imparted  to  the  object  by  the  blast  wave 
forces  would  tend  to  reduce  the  relative  velocity  between  the  object  and  the 
surrounding  gas.  Rotary  motion  imparted  to  the  object,  however,  would 
generally  result  in  an  increased  transfer  rate  at  given  locations  on  the  object. 
Whether  it  is  appropriate  to  consider  these  factors  in  greater  detail  depends 
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again  on  the  details  of  the  particulcr  application,  and  such  factors  are  not 
discussed  further  here. 

It  can  be  seen  from  Eq.  (6.3)  that  the  heating  durations  of  Figs.  6-1  and 
6-2  (of  either  the  bounding  or  reconunended  curves)  Increase  with  the  cube  root 
of  propellant  weight.  This  is  an  assumption  that  deserves  some  consideration. 
For  comparatively  small  propellant  quantities,  say  1000  lb  or  less,  the  fire¬ 
ball  duration  is  insufficient  for  appreciable  motion  (rise)  of  the  fireball, 
and  the  fireball  duration  is  then  essentially  synonymous  with  the  heating  dura¬ 
tion  oi  an  object  that  is  fixed  In  space.  For  larger  propellant  quantities, 
say  for  25,000  lb  and  more,  significant  motion  does  occur  and  the  heating  dura¬ 
tion  of  a  fixed  object  is  therefore  less  than  the  fireball  duration.  Thus,  the 
ratio  of  the  heating  duration  of  a  fixed  object  to  the  total  fireball  duration 
is  some  function  of  the  propellant  weight.  The  curves  of  Figs.  6-1  and  6-2  are 
based  on  measurements  fixed  in  space  at  the  25,000-lb  level,  and  extrapolation 
to  other  propellant  weight  levels  through  Eq.  (6.3)  inherently  assumes  an  invar¬ 
iance  of  this  ratio  of  durations.  For  application  to  weights  in  excess  of 
25,000  lb,  it  is  nevertheless  recommended  that  Eq.  (6.3)  be  used  in  conjunction 
with  the  curves  of  Figs.  6-1  and  6-2,  although  it  is  expected  that  the  curves 
would  be  somewhat  conservative.  For  extrapolation  to  significantly  lesser 
weights,  should  be  larger  than  given  by  Eq.  (6.3);  more  specifically,  at  the 
1000-lb  (or  less)  level,  as  given  by  Eq.  (6.3)  should  be  increased  by  a  mul¬ 
tiplying  factor  of  approximately  1.2  and  1.6  for  LO  /RP-1  and  LO  /LH  ,  respec- 
tively . 

It  is  possible  that  the  heat  transfer  hazard  can  be  intensified  by  the  oc¬ 
currence  of  chemical  activity  between  the  fireball  constituents  —  notably  the 
oxidants  —  and  the  surface  of  an  object  Immersed  in  the  fireball.  Predictions 
of  the  rates  (or  existence)  of  the  associated  chemical  reactions  are  not  inclu¬ 
ded  in  this  report,  in  part  due  to  the  heavy  dependence  of  such  reactions  on 
the  particular  application,  that  is,  on  the  molecular  constituents  of  the  ob¬ 
ject  and  the  surface  temperature  attained.  The  latter,  in  turn,  depends  on  the 
conf ig;uration  and  thermal  properties  of  the  object.  (The  reaction  also  depends 
critically,  of  course,  on  the  concentrations  of  various  atomic  and  molecular 
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species  —  and  their  excited  and  ionized  states  —  present  in  the  fireball.) 
Chemical  activity  is  mentioned  and  should  be  considered  in  any  application  —  par¬ 
ticularly  when  comparatively  lurge  propellant  quantities  are  involved  —  because 
tlie  reactions  can  provide  an  energy  contribution  (not  Included  in  Figs.  6-1  and 
6-2)  to  the  object. 

The  heat  flux  density  measuioments  upon  which  the  curves  of  Figs.  6-'.  and 
6-2  are  based  were  obtained  at  locations  no  closer  to  the  "center  of  explosion" 
than  about  one-fifth  of  the  radius  of  the  fireball,  and  it  would  be  expected 
that  the  heat  transfer  rates,  at  least  during  the  initial  "small"  fraction  of 
the  fireball  duration,  could  be  somewhat  more  severe  at  cr  "very  near"  the  cen¬ 
ter  of  explosion.  Passive  sensors  capable  of  providing  crude  indications  of 
comparatively  severe  heat  transfer  were  deployed  in  the  central  region  (within 
a  few  feet  of  the  planned  ignition  point)  throughout  most  ol  the  eleven  25,000- 
lb  tests  mentioned  above,  and  a  single  positive  Indication  was  obtained.  Spe¬ 
cifically,  from  0.1  to  0.2  in.  was  ablated  from  the  surface  of  a  solid  aluminum 

structure  in  such  a  way  as  to  suggest  comparatively  large  heat  flux  densities 

2 

over  limited  times,  for  instance,  of  the  order  of  1000  watt/ cm  for  2  sec.  (A 
thorough  analytic  evaluation  of  the  possible  ranges  of  heat  transfer  parameters 
resulting  in  the  above  ablation  has  not  been  performed;  for  details  ol  the  alum¬ 
inum  structure  and  its  ablation,  see  Appendix  C  of  Volume  1.)  It  is  not  clear 
if  chemical  activity,  as  mentioned  in  the  previous  paragraph,  was  an  energy  con¬ 
tributor  . 
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Appendix  A 

LIQUID  PROPELLANT  EXPLOSION  PHENOMENA  IN  RELATION  TO  TNT 

Tlie  characteristics  of  blast  waves  produced  by  explosive  energy  releases 
in  air  are  in  general  dependent  on  three  major  groups  of  properties:  those 
having  to  do  with  the  intrinsic  characteristics  of  the  explosive  material  (ex¬ 
plosive  properties);  those  having  to  do  with  the  manner  in  which  the  explosive 
material  is  assembled  (charge  properties);  and  those  having  to  do  with  the  en¬ 
vironment  surrounding  the  charge  (environmental  properties). 

Explosive  properties  would  Include,  for  example,  those  necessary  to  describe 
the  total  energy  release  per  gram  of  material  (that  contributes  to  the  explosion) 
and  the  rate  of  energy  release. 

Charge  properties  would  Include  those  necessary  to  describe  the  total  quan¬ 
tity  of  material,  its  shape,  and  any  confinement  effects  due  to  inert  material 
immediately  adjacent  to  the  charge. 

Environmental  properties,  for  example,  would  Include  those  necessary  to 
describe  the  nature  of  any  ground  surface  in  the  vicinity  of  the  charge,  the 
location  of  the  charge  with  respect  to  the  ground  surface  (height  of  burst), 
and  the  ambient  air  pressure. 

Ihe  explosive  properties  for  most  conventional  solid  high-explosive  mater¬ 
ials  are  quite  similar  and  for  a  given  material  are  often  adequately  described 
by  simply  stating  the  type  of  explosive.  For  liquid  propellants,  however,  the 
range  of  explosive  properties  and  thus  the  explosive  behavior  for  a  given  pro¬ 
pellant  combination  is  extremely  wide  and  generally  differs  considerably  from 
that  for  conventional  high  explosives. 

Specifically,  for  liquid  propellant  explosions  the  explosive  yield  or  frac¬ 
tion  of  the  total  energy  of  the  propellant  mixture  contributing  to  the  explosion 
is  not  only  dependent  on  propellant  type,  but  also  on  certain  other  parameters, 
including,  for  example, 
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1.  propellant  interaction  geometry  (mi':lng  mode) 

2.  time  of  Ignition 

3.  various  Initial  conditions,  such  as  spatial  and  velocity  distribution 
of  propellants  at  time  of  contact 


4. 

propellant  weight 

In 

distance 

addition,  the  shapes  of  the  propellant  pressure— distance 

curves  are  not  necessarily  identical  to  those  for  TOT. 

and  impulse— 

To  summarize,  the  parameters  which  generally  control  the  important  blast 

characteristics  of  liquid  propellant  explosions,  l.e.,  peak  overpressure  and 

♦ 

positive-phase  impulse,  are; 

1. 

total  propellant  weight 

W 

2. 

fraction  of  total  energy  contributing  to  explosion 

£ 

3. 

nature  of  reaction 

m 

4. 

ground  distance 

D 

5. 

height  of  burst 

h 

6  . 

charge  shape 

it 

By  means  of  normal  explosive  scaling  laws,  the  parameter  list  can  be  re¬ 
duced  to; 


1. 

nature 

of  reaction 

m 

2. 

scaled 

ground  distance 

D/(fW)^^^ 

3, 

scaled 

height  of  burst 

1/3 

V(fW)^' 

4. 

charge 

shape 

s 

For  some  conditions  it  is  possible  that  f  may  be  a  function  of  W. 


Confinement  effects,  ambient  overpressure,  and  nature  of  ground  surface  may 
have  to  be  considered  in  special  circumstances,  but  in  most  appllcatlms 
they  do  not  vary  sufficiently  to  require  consideration. 
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Considering  for  the  moment  a  surface  burst  condition  and  a  constant  charge 
shape,  only  two  basic  parameters  are  left; 

1.  nature  of  mixture  m 

1/3 

2.  scaled  ground  distance  D/(fW) 

Conceptually,  therefore,  all  desired  data  on  overpressure  (or  impulse)  for 
a  given  propellant  combination  (under  surface  burst  conditions  and  with  a  con¬ 
stant  charge  shape),  could  be  expressed  In  a  single  plot  as  follows: 


At  this  stage  It  is  appropriate  to  consider  in  more  detail  what  the  fac¬ 
tors  f  and  m  mean.  Earlier  f  was  defined  as  the  fraction  of  the  total  energy 
contributing  to  the  explosion.  This  actually  is  a  somewhat  loose  definition. 

It  could  be  argued,  for  example,  that  f  is  generally  not  the  effective  fraction 
of  energy  released  to  produce  a  given  overpressure  or  impulse  value  since  this 
also  depends  on  m,  as  illustrated  in  the  above  figure.  Ratt<er,  f  is  intended 
as  a  measure  of  the  total  amount  of  energy  released  in  the  explosion  that  ul¬ 
timately  contributes  to  the  blast  wave,  while  m  is  intended  as  a  measure  of  the 
manner  or  rate  of  energy  release.  The  reason  for  separating  the  factors  f  and  m, 


Considered  in  this  report  to  be  a  hemispherical  charge  on  the  surface. 
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rather  than  attempting  to  use  them  In  some  combined  fashion,  Is  that,  first, 
they  depend  on  somewhat  dlffei'ent  sets  of  basic  parameters  and,  secondly,  there 
is  generally  a  significant  range  of  distances  of  Interest  (at  intermediate  to 
long  distances  from  the  explosion)  where  the  m  values  become  unimportant  and  llie 
f  value  alone  controls  the  overpressure  and  Impulse  values.  Also,  in  this  re- 
gion  the  f  values  for  both  overpressure  and  Impulse  are  the  same. 

Within  the  definition  of  the  f  and  m  factors  given  above  there  is  still 
some  choice  of  how  to  quantitatively  determine  them  from  experimental  data  and 
to  use  them  in  a  prediction  method.  The  method  used  in  this  report  for  deter¬ 
mining  the  f  factor  is  based  on  TNT  equivalence  and  is  described  below. 

Assume,  for  example,  in  the  figure  below  that  curve  A  is  the  "standard"  TNT 
surface  burst  overpressure— distance  curve  and  curve  B  is  the  curve  for  a  par¬ 
ticular  teat  condition  involving  Wp  lb  of  propellant  (curve  B  points  are  plot¬ 
ted  using  Wp  for  computing  scaled  distance  values). 


This  occurs  because  as  a  blast  wave  proceeds  into  the  ambient  air  surround¬ 
ing  its  source,  its  behavior  tends  to  he  controlled  by  the  ambient  air  it  is 
propagating  through,  ana  the  Influence  of  the  specific  nature  of  the  source 
becomes  progressively  less  Important  with  increasing  distance  from  the  source. 
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For  each  <.  xpe.vimcnlal  point  on  curve  B  (P  D  )  the  weight  of  TNT,  WfTNT)  , 

1^1  1 

necessary  to  produce  the  same  overpressure  at  the  distance  as  obtained  in 
the  particular  test  is  computed  by: 

1  /3 

1.  detcimining  the  value  of  the  scaled  dlstcnce  (Dj^/W(TNT)j  '  )  from 

curve  A  ft  •  the  overpressure 

2.  dividing  by  this  scaled  distance  and  cubing  the  i-esults,  l.e.. 


WdNT)^ 


D^/W(TNT) 


3 


A  similar  process  waald  be  carried  out  for  the  impulse  data. 

The  ratios  of  W('INT)^/Wp  are  then  plotted  as  a  function  of  distance  as 
shown  below; 


The  f  factor  is  then  the  particular  value  of  the  ratio  of  W(TNT)^/Wp  ap¬ 
proached  by  both  the  overpressure  and  Impulse  data  at  long  distances.  Thus  the 
f  factor  is  simply  the  terminal  explosive  yield. 
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By  means  of  the  f  factor  so  determined,  the  experimental  data  points  can 
be  veplotted  using  fWp  for  computing  scaled  distance  and  compared  with  the  ideai 
curve  A  as  shown  below.  (C  is  curve  B  so  replotted.) 


A 


0/W 


1/3 


The  remaining  differences  between  curves  A  and  C  !  due  to  the  m  factor. 

The  PYRO  data  clearly  show  the  general  trends  discussed  above.  The  differ¬ 
ences  between  the  propellant  and  INT  overpressure— distance  and  impulse— dis¬ 
tance  curves  were  largest  for  the  hypergolic  combination,  and  both  pressure  and 
impulse  correction  factors  were  derived  as  a  function  of  scaled  distance.  For 
the  cryogenic  propellants  the  effects  were  somewhat  smaller,  and  it  was  only 
possible  to  derive  impulse  correction  factors.  The  p^’opellant  pressure  data 
tended  to  be  below  the  TOT  curve,  but  the  spread  was  Ip.rge  enough  so  that  the 
TNT  curve  was  used  as  a  conservative  bound.  It  should  be  emphasized  that  such 

a  bound  is  believed  to  be  overly  conservative  for  distances  closer  than  about  2 
1/3 

scaled  feet  (ft/lb  ).  Unfortunately  Insufficient  data  were  obtained  in  this 
region  to  actually  define  the  curve. 
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Appendix  B 

PEAK  OVERPRESSURE  VS  DISTANCE  AND  POSITIVE-PHASE  IMPULSE 
VS  DISTANCE  REFERENCE  CURVES 
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